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Abstract 

Electric field effects in iron- and copper-based superconductors were studied by using the 
first-principles calculations based on the density functional theory (DFT). The research objects 
include iron-based superconductors (KFe2Se2, LaFeAsO, NdFeAsO, and BaFe2As2) and 
copper-based superconductors (YBazCu307, HgBazCazCu30g, TleBazCaCu2Og, and 
Bi2Sr2Ca2Cu3010). To describe the strong correlation effect of 3d-electrons or 4f-electrons, the 
GGA+U method was used. Some results were further verified by the HSE method. The densities 
of states (DOS) were given. The change of the charge densities under electric fields is presented to 
demonstrate the electric field effect. It is found that the electron clouds of Fe ions in iron-based 
superconductors, Nd ions in Nd2Fe2As2O2, and Cu ions in copper-based superconductors change 
obviously. The pattern of the change is more like a rigid-body rotation than an elastic deformation. 
The author proposed that the rotation of the electron clouds of transition metal ions may be a new 
medium of superconducting electron pairing. The author’s views about some issues and 


suggestions on follow-up studies are also presented. 
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1. Introduction 

In 1986, J. G Bednorz and K. A. Miiller discovered that there may be high-temperature 
superconductivity (HTS) in La-Ba-Cu-O oxides"!, which revealed a new chapter in the research 
of superconductivity. In 1987, M. K. Wu et al. P! and Z. X. Zhao et al. ©! synthesized Y-Ba-Cu-O 
superconductors, respectively, almost at the same time. The critical temperature (Te) for 
superconductivity reached 90K and achieved a breakthrough above the liquid nitrogen 
temperature. In 1988, H. Maeda, Y. Tanaka, M. Fukutumi, and T. Asano. found BiSrCaCuO 
superconductors with a Ts above 100K “!. In 1988, Z. Z. Sheng et al. synthesized Tl-containing 
superconductors M6, In 1993, HgBazCuOu4zs, the first member of the Hg series, was synthesized 


by S. N. Putilin et al, with a Te of 94K [7]. In 1993, A. Schilling et al. synthesized 
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HgBa2CarCuzOs:s, raising the T, above 130K "!, The T, of Hg-1223 can reach 153 K under a high 
pressure Pl, In addition to copper-based superconductors, iron-based superconductors have also 
received widespread attention [101111112] The T, of PrFeAsOosoFo.11 reaches 52K [31] and 
SmFeAsOo.gs 55K l“, In 2012, Q. K. Xue et al. successfully grow FeSe thin films with one 
unit-cell thickness on the surface of SrTiO; substrates by the MBE method !!5!, The monolayer 
FeSe showed a T, above the liquid nitrogen temperature (77K). In 2014, Q. K. Xue et al. reported 
the FeSe monolayer with a Ts over 100K "41. 

The discovery of copper- and iron-based superconductors has seriously challenged the BCS 
theory, because BCS theory predicts that Te cannot be so high. Traditional electron-phonon 
interaction mechanism failed to explain the electron pairing mechanism and other new phenomena 
exhibited in copper- and iron-based superconductors. In 1987, P. W. Anderson l!” has put forward 
the famous RVB model for copper-based superconductors. Though many theories "897°! have 
been proposed so far, the electron pairing mechanism is still under intense debate. Anderson PH 
made it clear in 2007 that many theories about high-temperature superconducting electron pairing 
mechanism might be in the wrong direction. 

The author believes, to explain the electron pairing mechanism in high-temperature 
superconductors, we must go back to the traditional electron pairing mechanism. As early as 1950, 
H. Fröhlich 7! proposed that electron-lattice interaction is the source of superconductivity, 
supported by the isotope effect 7!?4!, Based on this, the well-known BCS theory ?*!?°! was 
developed, explaining the traditional superconductivity successfully. The author argues that the 
electronic-phonon interaction in BCS theory is closer to a mathematical tool than a simple and 
intuitive physical model. The simple and intuitive model is as follows. The electric field of an 
electron will distort the ions’ positions. When the electron leaves, the ions will not restore to the 
original positions immediately. This results in a positive charge region which will attract another 
electron and thus creates an attractive interaction between electrons. Based on the above 
consideration, the author believes that a study of the influence of electric fields on coper- and 
iron-based superconductors may provide important information about the electron pairing 
mechanism. Although many experiments of HTS under electric fields have been conducted 
27178129] the main purpose is to change the carrier concentration and thus the superconducting 
properties. It is related to new superconducting device preparation, no electron pairing mechanism 
involved. In the study of the electric field effect of eight typical unconventional superconductors, 
the author found, for the first time, the electron cloud of transition metal ions behaves like a 
rigid-body rotation. Based on this discovery, the author proposes a possible new superconducting 
electron pairing mechanism. This article describes the research methods, results, possible new 


superconducting electron pairing mechanism, and the authors' views on several issues about HTS. 


2. Research objects and calculation methods 

The author selected four typical iron-based superconductors (KFe2Se2, LaFeAsO, NdFeAsO 
and BaFe2As2) and four typical copper-based superconductors (YBazCu307, HgBazCa2Cu30s, 
Tl2Ba2CaCu20; and BizSr2CazCu30 10). It is worth noting that the objects calculated are not strictly 
actual superconductors. They are parent materials, and the author makes them closer to the actual 
superconductors by changing the stoichiometries or valence electron numbers. 

The first-principles calculation based on DFT was used. Most of the calculations were 
conducted within the Vienna Ab-initio Simulation Package (VASP) [3°] [31], The projector 
augmented plane wave (PAW) method !7!/°3] was used to describe the interaction between the 
valence electrons and the ion cores. The generalized gradient approximation (GGA) of 
Perdew-Burke-Ernzerhof (PBE) 4! was used for the exchange correlation potential. To describe 
the strong correlation of the 3d or 4f electrons, the GGA+U method is used. For the 3d orbitals of 
Fe and Cu, the Uer values are 4eV?>! and 6.5eV?*l7! | respectively. Plane waves were used to 
expand the valence electron wave functions. The cutoff energy is set to 400 eV which has been 
tested by taking both accuracy and computing cost into account. The valence electron 
configurations of the pseudopotentials are taken as 3p°3d’4s! for Fe (with a radius of 2.200 a.u.), 
4s*4p* for Se, 3s73p°4s! for K, 5s*5p°5d!6s? for La, 4s74p° for As, 4f'5s”5p°6s” for Nd, 5s*5p°6s” 
for Ba, 5d!°6s* for Hg, 3s?3pć4s? for Ca, 3p°3d!°4s! for Cu (with a radius of 2.000 a.u.), 2s?2p* for 
O, 4s74p°4d!5s? for Y, 5d!°6s’6p! for T1, 5d!°6s’6p* for Bi, and 4s74p°5s” for Sr. For some results, 
further calculations with the Heyd, Scuseria, and Ernzerhof (HSE) P3! hybrid-functional were 
performed for cross validation. 

Quantum Espresso (QE) simulation package P’! and norm-conserving pseudopotentials were 
employed in the calculation of iron-based superconductors under uniform electric fields. The 
valence electron configurations of the pseudopotentials are taken as 3s73p°3d°4s* for Fe (with a 
radius of 1.909 a.u.), 4s*4p* for Se, 3s73p°4s! for K, 5s*5p°5d'6s for La, 4874p? for As, 3s73p°4s! 
for Ba, and 2s2p* for O. Pseudopotentials with less valence electrons (3d°4s” for Fe, 4s! for K, 
and 6s? for Ba) were also used for comparison. The cut-off energy of the plane wave is set to 80 
Ry. The exchange-correlation interaction of electrons is escribed by the GGA+U method. Uerr is 
6.5 eV for the 3d orbital of Fe. Set occupations = 'fixed' in the input file. In general, occupations = 
'fixed' cannot be used for metal. Considering that the calculated objects are not good conductors, it 
is acceptable in the calculation of uniform electric fields, but it is difficult to converge to a high 
accuracy. In the following, the author will give an explanation. 

The results given below without description are obtained with VASP. When the QE software 
package is used, the author will point out. First, the geometry of the system is optimized until all 


forces on ions are less than 0.01 eV/A. After optimization, the influence of an electric field on the 


charge density is calculated. Calculations of iron- and copper-based superconductors were 
converged to 10° eV and 10°, respectively, for the tolerance on the total energy. There are three 
ways of applying an electric field to the system. The first is to insert one Li’ ion at an appropriate 
position (denoted by Li* Insertion); the second is to slightly change the position of some atom 
(Slight Change); and the third is to apply a uniform electric field (Uniform E-field) with QE 
package. 


3. Results and discussion 
3.1 K2Fe4Se4-2e and K2,Fe,Se4 

KFeSe superconductors have different stoichiometries, such as Ko.gFe1.6Se2!*°! and K,Fe2Se2 
M] The Te of KosFe16Se2 is 31.8K. To access an antiferromagnetic state, the author adopts 
KFe4Se4 for calculation. Two cases were calculated. One is K2Fe4Seq-2e, with two valence 
electrons subtracted, in which the valence of K is +1, Se is -2, and Fe is +2. The other is K2Fe4Seu, 
no electron added or subtracted, in which the valence states are rather complicated. If K is +1 and 
Se is -2, then Fe is +3/2. The k-point is set to 10x10x4 when using the GGA+U method, and 
6x6x2 when using the HSE method. Fig. 1 shows the DOS of K2Fe4Se,-2e, in which the 3d orbital 
of the Fe ion (Fe_d) is also presented. It’s an insulator, and the GGA+U method is in very good 
agreement with the HSE method except the bandgap. The GGA+U method gives a bandgap of 
about 1.3 eV, whereas the HSE method gives a bandgap of about 1.8 eV. 
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Fig. 1 DOS of K2Fe4Se4-2e from GGA+U (a) and HSE (b) method 
All the three methods were employed to study the influence of electric fields. The fractional 
coordinate of the inserted Li* ion is (0.5, 0.5, 0.62) and it is 2.55 A away from the nearest 
neighbor Fe, 2.73 A from Se, and 1.5 A from K. The original position the Se atom to be changed 
is (0.5, 0.5, 0.8682) and then revised to (0.5,0.5,0.8602). The position was changed in the (001) 
direction by -0.1 A. When applying the uniform electric field, QE package is used and the k-point 
is set to 6X62. The uniform electric field is along (111) direction, with a value of 0.00087 Ry a.u.. 


Set efield_cart(1) = 0.0005 Ry a.u., efield_cart(2) = 0.0005, and efield_cart(3) = 0.0005 in the 
4 


main input file. The convergence accuracy reaches 3x10 Ha. When the k-point is reduced to 
3x3x1, the convergence can reach 1x10% Ha, and the result is in accordance with that of 6x6x2 
k-point. 

Fig. 2 displays the crystal structure, spin density map and the charge density difference 
(CDD) with and without electric fields. The yellow color represents a positive value or increase of 
the charge density, while the blue color negative or decrease. It can be seen from the spin density 
map that there is a clear spin distribution on the iron ions and the antiferromagnetic state is present. 
For the three ways of exerting influence of electric field, the charge density around Fe ions 
changes significantly. There is no similar change around other atoms. The change is more like a 
rigid rotation rather than an elastic deformation, because the change is not entirely along the 
direction of the electric field. Some areas increase, while some areas decrease. Furthermore, the 
pattern of the change is like 3d electron clouds. The isosurface of the CDD is equal to that of the 
spin density, indicating the rotation is obvious. The author is sure that the observed results were 


not caused by low computational accuracy, from calculations with different accuracies. 
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Fig. 2 (a) Crystal structure of K2Fe4Se4 with one Li‘ ion inserted; (b) Spin density map (GGA+U, iso = 0.02 
num_electron/bohr’3); (c) CDD from Li* Insertion (GGA+U, iso = 0.02) ; (d) CDD from Li* Insertion (HSE, iso = 
0.02); (e) CDD from Slight Change (GGA+U, iso = 0.02); (f) CDD from Uniform E-field (QE, GGA+U, iso = 
0.02). iso represents isosurface. 

Table 1 shows the local magnetic moments (LMM) of the Fe ions in K2Fe4Se4-2e. When 
there is no electric field, LMM from the GGA+U method is very close to that from the HSE 
method, showing the consistency of the two methods. For the results from QE package, LMMs 
from the low accuracy and high k-point calculation are consistent with those from the high 
accuracy and low k-point calculation, indicating a converge of 3x10“ Ha is reasonable. The 
change of LMMs with and without electric fields is small, which indicates the electric field did not 


change the nature of the system significantly. 


Table 1 LMM (uB) of Fe atoms in K2FesSe4-2e 


Method Fel Fe2 Fe3 Fe4 
GGA+U No field 3.384 3.384 -3.384 -3.384 
Li* Insertion 3.453 3.322 -3.599 -3.274 
Slight Change 3.388 3.361 -3.385 -3.441 
HSE No field 3.376 3.376 -3.376 -3.376 
Li* Insertion 3.481 3.273 -3.507 -3.237 
QE(662k) No field 3.229 3.229 -3.193 -3.193 
Uniform E-field 3.276 3.301 -3.295 -3.299 
QE(331k) No field 3.227 3.227 -3.227 -3.227 
Uniform E-field 3.289 3.289 -3.289 -3.289 


In addition, increasing or decreasing the lattice constant by 2.5%, the same rotation can be 
observed. For Figure 2 (c), increasing the plane wave cutoff energy to 500eV, the results are 
consistent with that of 400eV. Calculations with valence electrons as 3d°4s7(1.9915 a.u.) for Fe 
and 4s! for K were performed. The k-point is 6x62 and the convergence reaches 1x10% Ha. The 
results are consistent with the previous, showing the observed rotation comes from the 3d 
electrons mainly. 

For K2Fe4Se4, with no electron subtracted or added, the results show that there is also a 
significant rotation of the electron cloud of the Fe ions. 
3.2 LaFeAsO 

The author calculated (LaFeAsO), directly, without increase or decrease of the number of the 
valence electrons. For the GGA+U method, the k-point is set to 10x10x6. For the HSE method, 
the k-point is set to 6x62. Fig. 3 shows the DOS. The two methods give similar DOS for the 


d-orbitals of the Fe atoms. 
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Fig. 3 DOS of (LaFeAsO)2 from GGA+U (a) and HSE (b) method 


When calculating the influence of electric field, the coordinate of the inserted Li* ion is (0.58, 


0.58, 0.62). The Li* ion is 1.877A, 1.979A, and 2.407A from the nearest neighbor Fe atom, As 
atom, and La atom, respectively. The original coordinate of the As atom to be changed is (0.75, 
0.75, 0.3363), and then revised to (0.75, 0.75, 0.3113). The position was changed in the (001) 
direction by -0.22 A. When using the QE package, the k-point is set to 6x6x2. The convergence 
reaches 1x105 Ha. Fig. 4 shows the crystal structure, CDDs and the total charge density (3d°4s” 


are taken as valence electrons for Fe). 
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Fig. 4 (a) Crystal structure of (LaFeAsO)2 with one Li* ion inserted; (b) CDD from Li” Insertion (GGA+U, iso = 
0.02) ;(c) CDD from Li‘ insertion (HSE, iso = 0.02); (d) Total charge density (GGA+U, iso = 0.3 ); (e) CDD 
from Slight Change (HSE, iso = 0.02); (f) CDD from Uniform E-field (QE, iso = 0.02) 

It can be seen from all the CDDs that there is a significant rotation of the electron clouds of 
the Fe ions. From Fig. 4(d), the distribution of the charge density around the iron atom does not 
have a perfect spherical symmetry. The author believes that this is the main reason for the rotation 
under electric fields. 

LMMs are given in Table 2. The results given by different methods are consistent, and the 
applied electric field has little effect on the local magnetic moments. 


Table 2 LMM of Fe in (LaFeAsO)2 


Method Fel Fe2 
GGA+U No field 3.261 -3.261 
Li* Insertion 3.286 -3.287 
Slight Change 3.280 -3.332 
HSE No field 3.282 -3.282 
Li* Insertion 3.244 -3.244 
Slight Change 3.349 -3.360 
QE No field 3.232 -3.222 
Uniform E-field 3.219 -3.230 


The calculation with 3d°4s* as valence electrons for Fe shows that there is a significant 
rotation of the electron cloud. Increasing the cut off energy of the plane wave to 500 eV gives the 
same results as 400 eV. The author also observed obvious rotation when 0.2 electrons added. 

In addition, the author calculated a 6x21 supercell. The position of the inserted Li* ion is 
(0.79, 0.38, 0.58); the coordinate of the As atom to be changed was (0.2083, 0.6250, 0.6637) and 
decreased by 0.01 (0.082 A) in the y direction. Fig. 5 shows the CDDs. The most important 
discovery is significant rotation occurs far away from the position where the electric field is 
applied. This shows that electron clouds of Fe ions in iron-based superconductors are more 
sensitive to the electric field, which is different from that of copper-based superconductors given 
below. The electron clouds of some Fe ions change obviously, but some not, which may be caused 


by the interaction between the charge density changes of different Fe ions. 
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Fig. 5 CDD of the 6x2 supercell from Li” Insertion (GGA+U, iso = 0.01) (a) and Slight Change (GGA+U, iso = 
0.01) (b) 
3.3 Nd2Fe2As202 

For the investigation of Nd2Fe2As202, the k-point is set to 10x10x5. Since Nd is a 4f 
transition metal atom, it has strong spin-orbital coupling and strong correlation. The spin-orbital 
coupling effect is considered in the calculation and the strong correlation effect is described by the 
GGA+U method (Uef = 7.0 eV). The calculation of uniform electric field together with 
spin-orbital coupling may be too expensive, so the author gives only the result of inserting one Li” 
ion. 

The Li* ion is located at (0.5, 0.25, 0.317), which is 1.912 A, 1.927 A, and 2.281 A from the 
nearest Nd, Fe, and As, respectively. The author also gives the results when the Li” ion is closer to 
one Nd atom. The coordinate of the Li* ion is (0.4, 0.4, 0.2) and it is 1.54 A from the nearest Nd 
atom. Fig. 6 shows the CDDs. When the distance between the Li ion and the Nd atom is large, 
rotation of the electron cloud of the Nd ion cannot be observed. When the distance decreased, 
obvious change appears. This suggests that the 4felectron cloud of Nd may also play an important 


role in superconductivity. 
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Fig. 6 CDDs from Li* Insertion (GGA+U, iso = 0.02). Li* ion is 1.912 A (a) and 1.54 A (b) from the nearest Nd 


atom, respectively. 
3.4 BarxK,FeAs2 

The T, of Baj..K,Fe2As2 can reach 38K 7!, Ba¿Fe4As4 and BaKFe4As4 were employed in the 
calculation, corresponding to x = 0 and 0.5, respectively. For the valence state of Ba2Fe4As4, As is 
-3, Ba is +2, and Fe is +2. For BaKFe4As4, As is -3, Ba is +2, K is +1 and Fe is +9/4. The k-point 
is set to 10104 when using the GGA+U method and 6x62 when using the HSE method. Fig. 7 
shows the DOS of BazFesAs4. The two methods are in good agreement, except that the GGA+U 


method gives a lower bandgap. 
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Fig. 7 DOS of BazFesAss from GGA+U (a) and HSE (b) method 
Two methods were used to calculate the electric field effect. One is to insert one Li” ion. The 
coordinate is (0.5,0.5,0.69), and it is 2.470 A, 2.191 A, and 2.268 A away from the nearest 
neighbor Ba, As, and Fe, respectively. The other is to apply a uniform electric field by the QE 
package, in which k-point is set to 662. The value and direction of the electric field are the same 
as those given above. The convergence reaches 1x104 Ha. Fig. 8 shows the crystal structure and 
the results, in which the CDD of BaKAs,Fe, from Li” Insertion is also given. The rotation of the 


charge densities of Fe ions is obvious. 
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Fig. 8 (a) Crystal structure of Ba2Fe4As4; (b) CDD of BarFesAsa from Li” Insertion (GGA+U, iso = 0.02); (c) CDD 
of Ba2Fe4As4 from Uniform E-field (QE, GGA+U, iso = 0.02); (d) CDD of BaKAsaFea from Li* Insertion 
(GGA+U, iso = 0.02) 

For BaKAs4Fe4, apparent charge density rotation can also be observed when the position of 
an As atom is slightly changed or an uniform electric field is applied 
3.5 (YBa2Cu307)2 

To access an antiferromagnetic state, (YBazCu307)2 was used, no electron added or 
subtracted. When the GGA+U method is adopted, the k-point is set to 8x8x4. When the HSE 
method is used, the k-point is 4x4x2. Fig. 9 shows the DOS. 
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Fig. 9 DOS of (YBaxCus07)2 from GGA+U (a) and HSE (b) method. Cu_d has been multiplied by 3 for clarity. 

The spin density of the system was calculated. The effect of an electric field induced by a 
slight change in the position of one oxygen atom in the copper oxide plane was calculated. The 
original coordinate of the oxygen atom was (0.75, 0.25, 0.62245), and then changed to (0.8, 0.3, 
0.62245). The position was changed in the (110) direction by 0.39 A. The method is GGA+U, and 
the result was validated by HSE. The author also performed the calculation for the vertical change 
(0.23 A) of the oxygen atom. Substituting Cu with Zn was also calculated with the GGA+U 
method (3d electrons of Zn are regarded as valence electrons and set Uet = 6.5eV), and the 


position of oxygen atom was changed along (110). Fig. 10 shows the crystal structure and results. 
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(d) 
Fig. 10 (a) Crystal structure; (b) Spin density map (GGA+U, iso = 0.01); (c) CDD form Slight Change (GGA+U, 
iso = 0.004); (d) CDD from Slight Change (HSE, iso = 0.004); (e) CDD from Slight Change along z direction 
(GGA+U, iso = 0.004); (f) CDD of YBaZnO from Slight Change (GGA+U, Uef = 6.5eV, iso = 0.004) 

It can be seen from the spin density map that the Cu ions in the copper-oxygen layer exhibit 
an antiferromagnetic state. The electron cloud of Cu ions shows an obvious rotation. It is worth 
noting that the isourface is about one order of magnitude smaller than that of Fe ions in iron-based 
superconductors. The electron cloud of Zn did not rotate at all. Calculations with Uefe = 0 show 
that the electron cloud of Zn does not rotate too. Table 3 gives the LMM of the Cu atom. When the 
electric field is not applied, the results given by the GGA+U and the HSE method are very 
consistent. When there is an influence, the change of the LMM is very small. 


Table 3 LMM of the Cu atom in (YBa2Cu307)2 


Method Cu 
GGA+U_ No field 0.634 

Slight Change xy 0.631 

Slight Change z 0.641 
HSE No field 0.631 


The rotation is also observed in (YBazCu307)2 with 0.3 or 2 electrons added. The author 
believes that, if the 3d orbitals of Cu is not fully filled, the electron will rotate under electric fields. 
When the 3d shell is filled, the rotation cannot be observed. 

3.6 (HgBa2Ca2Cu30s)2 

To achieve an antiferromagnetic state, (HgBazCa2Cu30s)2 was used. No electrons were added 
or subtracted. With the GGA+U method, the k-point is set to 8x8x4. With the HSE method, the 
k-point is set to 4x42. Fig. 11 shows the DOS. The 3d orbital of Cu (Cu_d) is of the middle 
copper-oxide layer. 
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Fig. 11 DOS of (HgBa2Ca2Cu3Os)2 from GGA+U (a) and HSE (b) method. 

The effect of electric fields is calculated by slightly changing one oxygen atom in the middle 
copper-oxygen layer. The original position of the oxygen atom is (0.75, 0.25, 0.5) and then 
changed to (0.8, 0.3, 0.5). The position was changed in the (110) direction by 0.39 A. The result 
was verified by the HSE method. Fig. 12 shows the crystal structure, spin density map and CDDs 
with and without electric fields. The charge density of Cu ions exhibits clear rotations and the 


GGA+U method gives consistent result with the HSE method. 
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Fig. 12 (a) Crystal structure; (b) Spin density map (HSE, iso = 0.01); (c) CDD from Slight Change (GGA+U, iso = 


0.003); (d) CDD from Slight Change (HSE, iso = 0.003) 
Table 4 gives the LMM. Consistency between the HSE and GGA+U method can be seen and 
the effect of a Slight Change on the LMM is very small. 


Table 4 LMM of the Cu atom in (HgBa2Ca2Cu30s)2 


Method Cu 
GGA+U No field 0.603 
Slight Change 0.600 
HSE No field 0.618 
Slight Change 0.619 
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Calculations of (HgBazCa2Cu309)2/(HgBazCa2Cu309+2e)2 show that there are significant 
rotations of the electron cloud of copper atoms. 

Further, a 4x11 supercell of (HgBazCa2Cu30s)2 was calculated. The original coordinate of 
the oxygen atom to be changed is (0.1875, 0.25, 0.5), and then changed to (0.2, 0.3, 0.5). The 
CDD is displayed in Fig. 13. There is a clear rotation on Cu atom. However, no rotation occurs far 
away from where the influence is applied. This is totally different from that of iron-based 


superconductors. For iron-based superconductors, rotation can be observed further away. 


Fig. 13 CDD of the 4x11 supercell from Slight Change (GGA+U, iso = 0.003) 
3.7 (ThBazCaCu2Os)4 
To access an antiferromagnetic state, the author adopted (Tlz2BazCaCuz2Os)4, no electron added 
or subtracted. The GGA+U method is used, and the k-point is set to 6x6x1. Fig. 14 shows the 
DOS. 
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Fig. 14 DOS of (ThBaxCaCu20s)s from GGA+U method 
The position of the O atom was changed by 0.39 A, -0.39 A and 0.78 A, respectively, in the 
(110) direction. Fig. 15 displays the CDDs. It can be clearly seen that the electron cloud of the 
Cu? ion rotates with the change of the O atom’s position. Slightly changing the position of 
oxygen atoms near TI or Ba, the author did not observe the rotation of the electron cloud of TI or 


Ba. 


(a) (b) 


Fig. 15 (a) Copper-oxygen layer; (b) CDD from Slight Change of 0.39A; (c) CDD from Slight Change of -0.39 A; 
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(c) CDD from Slight Change of 0.78A (GGA+U, iso = 0.003) 
3.8 (BizSr2Ca2Cu3O10)4 
(Bi2Sr2Ca2Cu3010)4 was adopted, no electron added or subtracted. The k-point is set to 5x5x1. 
Fig. 16 shows the DOS. 
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Fig. 16 DOS of (Bi2SrxCa2Cu3O10)4 from GGA+U method 
The system has three adjacent copper-oxygen layers, the position of one oxygen atom in the 
uppermost copper-oxygen layer is changed. The original coordinate is (0.75, 0.25, 0.586189), and 
then changed to (0.8, 0.3, 0.586189). The author also studied the case of increasing or decreasing 
the lattice constant. Fig. 17 shows the results. There is also a significant rotation of the electron 


cloud of Cu ions when the lattice constant is increased or decreased. 


Fig. 17 (a) CDD from Slight Change (GGA+U, iso = 0.002), (b) The lattice constant decreases by 2.5%; (c) The 
lattice constant increases by 2.5% 

It can be seen from the results of the above eight systems all the electron clouds of Cu ions or 
Fe ions have obvious changes under electric fields.The author believes this change should be 
called a rotation not an elastic deformation, because the CDD does not simply follow the direction 
of the electric field, but increases in some areas and decreases in others. 

Is this rotation just of 3d (or 4f) electron cloud? Considering the orthogonality of the atomic 
orbitals, the author thinks a change of the outer electron clouds will inevitably affect the inner 
electron clouds. So, the inner electron cloud will rotate too. However, the rotation of the inner 
electron cloud is difficult to observe because they have spherical symmetry. 

Why do electron clouds rotate? We know that transition metal ions (3d or 4f) have 
incompletely filled 3d or 4f orbitals. The electron cloud does not have perfect spherical symmetry. 
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It is like a polar molecule and will rotate under electric fields. 

In addition to the above eight typical unconventional superconductors, other superconductors, 
such as metallic Nb, heavy fermion superconductor UPt:, rare earth superconductor SmRh4B,, et 
al., have incompletely filled 3d or 4f shells. They are likely to rotate under an electric field and 


thus play an important role in superconducting electron pairing. 


4. Possible Superconducting Electron Pairing Mechanism, Views on some issues, 
and Research Prospect 
4.1 Mechanism 

Mr. Anderson has raised the issue of whether a glue exists in high temperature 
superconductors. The author believes it does exist. Combined with the simple and intuitive model 
of electron-phonon interaction and the fact that transition-metal ions’ electron clouds rotate under 
an electric field, the author proposes that the rotation of the transition metal electron cloud is the 
glue. When an electron comes to a new location, the electron cloud of the transition metal ions in 
the vicinity will rotate. In this way, the charge density around the electron will decrease 
(equivalent to the emergence of a positive charge). When the electron leaves, the electron cloud of 
the transition metal ions will not relax immediately (it is so crucial that that the author will discuss 
it later), so that there will be a region lack of charge (equivalent to a positive charge region) that 
will attract another electron. Attraction between electrons appears. This mechanism is essentially 
the same as the electron-phonon interaction, except that the medium is the rotation of the electron 
cloud, not the displacement of ions. 
4.2 Author’s views on some issues 
4.2.1 Relaxation time of the electron cloud rotation 

The most important point in the new pairing mechanism is the relaxation time of the electron 
cloud's rotation. If the relaxation time is smaller than that of the elastic deformation of the electron 
cloud, then it cannot be the glue. The author believes that this relaxation time will not be so small, 
because the inner electrons will also rotate due to the requirement of orthogonality of electron 
wavefunctions. This is the most important follow-up study. 
4.2.2 Is the rotation of electron clouds of ions or ions? 

The author think ‘Rotation of ions’ is better, because it is more simple and intuitive. 
4.2.3 The role of strong correlation 

For the electron cloud rotation, iron-based superconductors are more obvious than 
copper-based superconductors, but the T. of copper-based superconductors is higher than that of 
iron-based superconductors. In addition to different electronic structures, the author believes that 


the strong correlation plays an important role. Copper-based superconductors have stronger 
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correlation effect than iron-based superconductors. The screening is strong in iron-based 
superconductors, which causes a lower T.. For the relationship between strong correlation and 
high Te, the author believes strong correlation is a necessary but not sufficient condition for high 
Te. In strongly correlated systems, electrons cannot be described by free electrons or weakly 
interacting quasi-particles, so that a strong attraction is possible. The screening in traditional 
superconductors is strong, on the one hand reducing the repulsion between electrons, but also lead 
to a weak attraction. 
4.2.4 Relationship between magnetism and unconventional superconductivity 

The author thinks that the rotation of the electron cloud of transition metal ions is the 
fundamental reason for the unconventional superconductivity. The rotation is because the atomic 
shells are not fully filled and no spherical symmetry present. When there is an electric field, it will 
rotate. Magnetism is also because the shells are not fully filled. 
4.3 Research prospect 
4.3.1 The research in this article is mainly concentrated on the existence of the electron cloud 
rotation. More detailed researches are needed. For example, the projected DOS and more analysis 
should be given. In addition, the study also need more accurate calculations, such as the use of the 
Quantum Monte Carlo method. Most results of this paper are qualitative. Further research requires 
quantitative analysis, such as quantification of the applied electric field and quantification of the 
electron cloud rotation. 
4.3.2 The author believes that determination of the relaxation time of the electron cloud rotation is 
crucial. Theoretical and experimental studies of the electron cloud rotation and the relaxation time 
should be given. 
4.4.3 New theories and application. For example, theory of the interaction between electrons and 
the rotation of ions, possible long-range modes of the rotation in iron-based superconductors, and 


influence (or application) of rotation of transition metal ions in other fields. 


4. Conclusion 

The author studied eight typical unconventional superconductors. The electron cloud of 
transition metal ions is found to rotate under electric fields. Extensive calculations show that if 
there is an incompletely filled electronic orbital, the electron cloud will rotate. If the orbital is 
filled, the electron cloud will not rotate. The authors believe that the rotation of transition metal 
ion’s electron cloud may be a new medium of superconducting electron pairing, just as lattice 
deformation is the conventional superconducting electron pairing medium. The author’ views on 
this new possible pairing mechanism, strong correlation and superconductivity, magnetism and 


superconductivity, and follow up studies are given. 
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